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A variational technique is established to deal with the Schrodinger equation with parity-time (VT) symmetric 
Gaussian complex potential. The method is extended to the linear and self- focusing and defocusing nonlinear 
cases. Some unusual properties in PT systems such as transverse power flow and PT breaking points can be 
analyzed by this method. Following numerical simulations, the analytical results are in good agreement with 
the numerical results. 
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INTRODUCTION: People have paid much attention 
to optical properties in parity-time {VT) symmetric 
system during recent years both theoretically and ex- 
perimentally [1-6]. Bender et al found that the non- 
Hcrmitian Hamiltonian with VT symmetry can exhibit 
entirely real spectrum [7]. In optics, VT symmetric po- 
tentials can be also constructed by a complex gain/loss 
refractive-index distribution into the waveguide, which 
make the refractive index distribution is VT symme- 
try [1-3]. Many unusual features stemming from VT 
symmetry have been found, such as power oscillation [8], 
absorption enhanced transmission [6], nonlinear switch- 
ing structure [9], and unidirectional invisibility [10]. 

But the properties of VT symmetry mentioned above 
arc mostly discussed by numerical simulations and ex- 
perimental results. In order to get a good understanding 
of the properties of VT symmetry solitons, it is essen- 
tial to present an analytic solution, even an approximate 
one. Musslimani et al have present closed form solutions 
to a certain class of nonlinear Schrodinger equations in- 
volving VT symmetry potential [11]. R.EI-Ganainy et al 
have developed a formalism for coupled optical parity- 
time symmetric systems by Lagranian principles [1]. It 
would be very desirable to obtain approximate analyt- 
ical results for the general VT symmetry Schrodinger 
equations. For the nonlinear Schrodinger equations, the 
approximate solutions can provide a better physical un- 
derstanding of the interplay between VT symmetry po- 
tential and nonlinear effects in connection with solitons 
propagation. The variational method has been widely 
applied to obtain approximated solutions for problems 
concerning beam propagation within the framework of 
the NLSE [13]. Motivated by the previous works in which 
it was studied the dissipative solitons solutions by the 



variational approach [14, 15]. In this paper, we will em- 
ploy a variational approach based on the dissipative sys- 
tem to deal with the Schrodinger equation with PT sym- 
metric Gaussian complex potential, and some unusual 
properties in VT system can be analyzed by this method. 
MODEL: We consider the (l+l)-dimensional evolu- 
tion equation along the longitudinal direction z in lin- 
ear or Kerr-nonlinear media with complex PT poten- 
tials. The propagation of the beam is governed by the 
Schrodinger equation [1], i.e. 

id:,U + dlU + [V{x) -t- iW{x)]U + ct\U\^U = 0. (1) 

Here U is the slowly varying complex field envelop, 
X is the transverse coordinate, a = 1 represents the 
self- focusing propagation, ct = — 1 represents the self- 
defocusing propagation, and a = Q represents the linear 
situation, ^(a;) and W{x) are the real and the imaginary 
parts of VT symmetry complex potentials respectively. 
The normalization relation is given in the reference [1] . 
If the solution of stationary wave modes for Eq. (1) has 
the form like U = (j){x) cxp[i6'(a;)] exp(z6z), where b is the 
propagation constant, 4>{x) and 0{x) are real functions 
standing for the amplitude and phase of light, then 

dlcf> + V{x)^ + a\(j>\^cf>-b^^cf>id^ef (2) 

d^{(j>^d^J) + W{x)^^ = 0. (3) 

It is noteworthy that Eq. (3) is a specialized form of 
Poynting's theorem for electromagnetic fields in disper- 
sive media with losses [12], where the transverse com- 
ponent of Poynting vector can be written as S* = 4)^3x8 
(only the transverse differential remains for stationary 
modes). The term W{x)(jP represents the absorption(or 



gain) in the medium. Equation (3) can be applicable to 
the general dissipative system. At the center of the VT 
system, we have 
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where a(z), w{z)^ C{z), and '^{z) stand for the ampli- 
tude, width, chirp and the phase of the beam, 9q{z)^ and 
Wi(z) stand for the amplitude and the width of wavefront 
lean 6. By substituting the trial solution into Eq.(5) and 
the averaged conservative Lagrangian can be obtained. 
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The Eulcr-Lagrangc equations for parameter rj is 
d f d < Lc >\ d < Lc > 
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in which Q = -'iW{x)U is the dissipative term in Eq. 
(1). Following the process [14,15], substituting rj with a, 
w, wi, 6q, C and $ into the Eq. (8), we arrive at the 
relations between the parameters. 
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(4) 
That means the energy-flow at the center must equal 
to the total loss (or gain) at each side. As a result, the 
intensity at the center can not be zero for any stable 
solution and the dark solitons can not exist in any VT 
system. 

VARIATIONAL METHOD: 

According to the literatures [14,15], the Lagrangian 
for Eq. (1) is the sum of two terms, a conservative one 
Lc and a non-conservative one L^c, where 

Lc = '-{Ud,U'-U'd,U) + \d.U\'~V\U\^-a^\U\\ (5) 

The VT symmetric gaussian complex potentials are 

used here as V{x) = Vbexp(— a;^) and W{x) ~ Ml'^j 

Woxexp(—x2), where Vq and Wo represent the depth of 

the real and the imaginary parts of VT potentials re- b 

spcctively [16]. We suppose the trial function as. 



Equations (9) is a trial result which means the power of a 
stationary modes is unchange during propagation. From 
Eqs.(lO) and (11), we find that the shape of the phase 
function 9 is independent of the beam width. Numeri- 
cal simulations show that even for unstable solution, the 
shape of the phase function is very similar to the error 
function, and its amplitude (^o) is mainly determined by 
Wo- We also find that these properties are identical to 
the multimodc solitons and other VT systems such as 
Scarff II potential. 

Introducing Pq = ^/tto^w, we arrive at the evolution 
equations. 
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LINEAR MODES: For the case of stable solution, the 
derivatives of beam width and curvature in the evolution 
equations are vanished. In linear propagation, i.e. a = 0, 
we can get (7 = 0, n{w) = 0, 
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From Eqs. (16), for a large value of Vq, we have a max- 
imum value of Woci which can be regarded as the VT- 
brcaking point. For a small value of Vq, there doesn't 
exist a maximum value of Woe, however, when Wq = 
2^/^y/Vo, 6-^0, and W^ can be regarded as the VT- 
breaking point. Figures (l)(a) and (b) shows the relation 
to the eigenvalues with Wo for different Vq in numeri- 
cal and variational method. We can see that for a large 
K) (1^=10), the Woe obtained by the variational method 
is in good agreement with the 7^7~-breaking point in the 
numerical method, and for a small Vo (Vq = 2), the W^ is 
close to the PT-breaking point. The PT-brcaking point 
with different Vo in numerical method and the relation 
of the Woe or the Wq with Vq in analysis are shown 
in the fig. (l)(d). One can see that the Woe or the W^ 
is in good agreement with the T^T-breaking point. Fig- 
ure(l)(c) shows that the amplitude and the nontrivial 
phase of the eigenfunction are in good agreement with 
the numerical results too. 

SOLITONS: For the nonlinear propagation, i.e. a = 
±1, we can get C = and aPo = 2-\/27rw'^/i(u;) for the 
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Fig. 1. (color online) (a) and (b)The eigenvalues for lin- 
ear VT modes versus the value of Wq for (a) Vo ~ 2 and 
(b) Vq = 10. (c) The amplitude and the nontrivial phase 
of the eigenfunetion with Vq = 10 and Wq ^ 10 {d)'PT- 
breaking point versus Vb-(For all of the cases, solid: vari- 
ational results, dotted: numerical results) 



VT solitons. For a given value of Vq and Wq, the VT 
solitons obey 
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Without losing of generality, we take Vb = 10 and Wo = 
5 in this subsection. 

Figures (2) (a) and (b) are the relations between the 
power and the beam width of solitons with the propa- 
gation constants, respectively. The critical propagation 
constant be at zero power is the eigenvalue for the linear 
mode. When propagation constant b > be, it is corre- 
sponding to the solitons in the self-focusing nonlinear 
media, otherwise it represents the solitons in the self- 
defocusing media. We can see that in the self-defocusing 
media, the power of the bright solitons decreases with 
increasing of the propagation constants, however, it is 
contrary to the bright solitons in the self-focusing me- 
dia. However, the beam width decreases with increasing 
of the propagation constants in the whole region. The 
amplitude and the nontrivial phase of the solitons in nu- 
merical method and variational approach arc shown in 
Figs. (2)(c) and (d). which are corresponding to the cases 
marked as circle symbols in Figs. 2(a) and (b). One can 
see that the analytical results are in excellent agreement 
with the numerical results for both self-focusing and self- 
defocusing cases. 

Stability of the solitons: From the evolution equations, 



Fig. 2. (Color online) (a) The power of the solitons versus 
propagation constants; (b) the beam width of the soli- 
tons versus propagation constants; (Red: self-defocusing 
nonlinear, blue: self- focusing nonlinear. ); (c) and (d) 
The amplitude and the nontrivial phase of the solitons in 
self-defocusing and self-focusing nonlinear , respectively, 
which is corresponding to the cases marked as circle sym- 
bols in (a) and (b). (Red: the nontrivial phase of the soli- 
tons , blue: the amplitude of the solitons. For all of the 
cases, solid: analytical results, dashed: numerical results, 
and Vb=10, Wo=5. ) 



Eq. (18), we define 
F 
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The equation associated with the the Jacoby determi- 
nant corresponding to the stability criterion is [14] 



A^+4CA-32C2-4z«^-^^ 
aw 



2ttw^ 



(19) 



When Re{X) = 0, or A^ < 0, solitons is stable. Substi- 
tuting C = and ctPq = 2y/2Trw^ fi{w) to Eq. (19), one 
has 

A2 ^ 4w-^^ + 12n(w) < 0. (20) 

aw 

CONCLUSION: In conclusion, we have analyzed the 
properties of the beam in FT symmetric Gaussian com- 
plex potentials by the variation method. According to 
the analysis, the transverse power-flow density is associ- 
ated with the wavefront lean 6 of the beam. For the linear 
case, the T'T-breaking points are obtained by the varia- 
tional approach, and the eigenvalue obtained by the vari- 
ational approach are in agreement with the ones by the 
numerical method. For the nonlinear case, there exists a 



critical propagation constant for bright soliton existing 
in the self-focusing and self-defocusing media. The rela- 
tions between the power of solitons with the propagation 
constants and the beam width are gotten analytically 
and confirmed by numerical simulations. The amplitude 
and the nontrivial phase of the solitons in both the self- 
focusing and self-defocusing media are agreement with 
the numerical results. 
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